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Abstract 

Conserving approximations are applied to the attractive Holstein and Hub- 
bard models (on an infinite-dimensional hypercubic lattice). All effects of 
nonconstant density of states and vertex corrections are taken into account in 
the weak-coupling regime. Infinite summation of certain classes of diagrams 
turns out to be a quantitatively less accurate approximation than truncation 
of the conserving approximations to a finite order, but the infinite summation 
approximations do show the correct qualitative behavior of generating a peak 

in the transition temperature as the interaction strength increases. 
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I. INTRODUCTION 



It is generally believed that the theoretical aspects of conventional superconductors are 
well understood and that quantitative predictions agree with experimentS'i. The reason 
why low temperature superconductors can be described accurately for all physical values of 
the electron-phonon coupling is due to Migdal's theorem!: the ratio of the electron mass 
to the ion-core mass provides the small parameter that guarantees rapid convergence of 
the theory. EliashbergS generalized Migdal's theorem to the superconducting state and 
provided the framework for quantitative calculations of the superconducting properties of 
real materials!!!. 

A more precise way of stating Migdal's theorem is to say that only the electrons that lie 
in an energy shell of width fioebye about the Fermi surface are affected by phonon scattering, 
and the only important scattering events involve the virtual emission and reabsorption of 
phonons in an ordered fashion, where the last emitted phonon is the first absorbed phonon, 
and so on. Migdal-Eliashberg (ME) theory neglects vertex corrections (which involve cross- 
ings of the phonon lines) and is an accurate approximation for small phonon frequencies. 
The remaining unanswered question is how large does the phonon frequency have to be 
before the effects of vertex corrections are observable? 

There are many materials that are hypothesized to be electron-phonon mediated su- 
perconductors, but have large phonon frequencies. Bai_ x K x Bi03 is a charge-density-wave 
(CDW) insulator at zero doping (x = 0), but becomes a superconductor (SC) away from 
half-filling! (x > 0.37). The maximum phonon frequency is fi max . = 80 meV for the optical 
oxygen modes,! while the bandwidthB is W = 4 eV, so the ratio of the phonon energy scale 
to the electronic energy scale is VL^-^JW = 0.02. Are vertex corrections important for this 
material? 

Alkali-metal-doped Ceo is another superconducting material that is hypothesized to have 
an electron-phonon pairing mechanism. There are very high frequency phonons that cor- 
respond to distorting the C 60 cage! (Q m 0.2 eV), while the electronic bandwidth is quite 
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narrowu (W = 1.0 eV), resulting in £l ma , x ./W ~ 0.2. Clearly, vertex corrections must play 
an important role in phonon-mediated pairing mechanisms in these materials. 

The effect of vertex corrections on superconducting properties, in particular, on the 
superconducting transition temperature, have been studied in the past0"0. Grabowski and 
Sham0 showed that vertex corrections lower T c for the repulsive electron gas with T c — > 
for some critical value of the plasma frequency. The electron-phonon interaction has also 
been examined,0Jii and, in general, vertex corrections also cause T c to drop. Is this always 
the effect of vertex corrections, or can vertex corrections sometimes cause an enhancement 
toT c ? 

In this contribution the effects of vertex corrections are examined in a systematic fashion 
via weak-coupling conserving approximations for the attractive HolsteinEl and Hubbard0 
models. The effects of Coulomb repulsion are explicitly neglected here. A detailed com- 
parison of these perturbation schemes can be made to exact results for these models in 
the infinite-dimensional limit to determine which weak-coupling approximation is the most 
accurate. 

The Holstein model consists of conduction electrons that interact with localized (Ein- 
stein) phonons: 
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H = -7Tn E ( c Uka + clc Ja ) + 5>x, - n)(nn + n n - 1) + ^Mfi 2 £ x) + ~ £ ^ 

2 ^ d (j,k)a j Z j 1 j M 

(1) 

where cJ CT (cj a ) creates (destroys) an electron at site j with spin a, rij a = c^ a Cj a is the electron 
number operator, and Xj (pj) is the phonon coordinate (momentum) at site j. The hopping 
matrix elements connect the nearest neighbors of a hypercubic lattice in (i-dimensions. The 
unit of energy is chosen to be the rescaled matrix element t*. The phonon has a mass M 
(chosen to be M = 1), a frequency Q, and a spring constant k = MQ 2 associated with it. 
The electron-phonon coupling constant (deformation potential) is denoted by g; the effective 
electron-electron interaction strength is then the bipolaron binding energy 
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The chemical potential is denoted by fi and particle-hole symmetry occurs for \i = 0. 

In the instantaneous limit where U remains finite and g and Q are large compared to the 
bandwidth (g,Q — ► oo,U = finite), the Holstein model maps onto the attractive Hubbard 
model0 

t* 11 

H = -TTn E ( c 3a + c Lc jCT ) -/iE(% + n n) + u XX n iT - 7;)( n 3i ~ ~) ( 3 ) 

with U defined by Eq. (@). 

The weak-coupling theory is based upon the conserving approximations of Baym and 
Kadanofifi the free energy functional $ is approximated by a series expansion of skeleton 
diagrams of the dressed Green's function G; the self energy £(zu; n ) is determined by func- 
tional differentiation E(iu; n ) = 5<& / 5G(iu n ) at each Matsubara frequency uo n = (2n + l)vrT; 
and the irreducible vertex functions Y(iuj mi iuo n ) (in the relevant channels) are determined 
by a second functional differentiation. 

Independently, Van Donge n§ and Martin-Rodero and Floreslll showed that the free 
energy must be expanded to order U 2 to determine the correct transition temperature 
in the limit \U\ — > for the Hubbard model at half-filling. The vertex corrections re- 
duced the Hartree-Fock transition temperature by a factor of order three, but the gap ratio 
2A(0) /ksT c fa 3.53 was unchanged (to lowest order). 

In addition to reproducing the weak-coupling limit properly, one hopes that the conserv- 
ing approximations will also be able reproduce the peak in the transition temperature as a 
function of interaction strength that occurs as the system crosses over from a weak-coupling 
regime (where pair formation and condensation both occur at T c ) to a strong-coupling regime 
(where preformed pairs order at a lower temperature It will turn out that this feature 
is not easily reproduced by a truncated conserving approximation. 

The infinite-dimensional limit of Metzner and Vollhardt0 is taken (d — > oo), in which 
the electronic many-body problem becomes a local (impurity) problem that retains its com- 



plicated dynamics in time. The large-dimensional limit is quite useful because both the 
Holsteini and Hubbard&H models can be solved exactly using the quantum Monte Carlo 
(QMC) techniques of Hirsch and These exact solutions have many of the qualitative 

features of the many-body problem in finite dimensions. They also provide a unique testing 
ground for various weak-coupling theories, since the approximate theory can be compared 
directly to the exact solution in the thermodynamic limit. 

In the infinite dimensional limit, the hopping integral is scaled to zero in such a fashion 
that the free-electron kinetic energy remains finite while the self energy for the single-particle 
Green's function and the irreducible vertex functions have no momentum dependence and 
are functionals of the local Green's functionl^'BEl. This limit retains the strong-correlation 
effects that arise from trying to simultaneously minimize both the kinetic energy and the 
potential energy. 

The many-body problem is solved by mapping it onto an auxiliary impurity problemllil 
in a time- dependent field that mimics the hopping of an electron onto a site at time r and 
off the site at a time r' . The action for the impurity problem is found by integrating out all 
of the degrees of freedom of the other lattice sites in a path-integral formalism.^ The result 
is an effective action 

Se ff . = EaI P dTj p dT'cUT)GoHr ~ r')c a (r') + E<r So dr[gx{r) - ^M(r)c a (r) - 1] 

+ \M^dT[n 2 x 2 {T)+x 2 {r)} (4) 

where Gq 1 is the "bare" Green's function that contains all of the dynamical information of 
the other sites of the lattice. The interacting Green's function, defined to be 

is determined by Dyson's equation 

G- 1 = G-\iu n ) = Go\iu n ) - £(iw n ). (6) 

A self-consistency relation is required in order to determine the bare Green's function 
Gq. This is achieved by mapping the impurity problem onto the infinite-dimensional lattice 



thereby equating the full Green's function for the impurity problem with the local Green's 
function for the lattice 

G j:j (iuj n ) = G(k, iuj n ) = Y^i iuJ n + /i - E(k) - E^o^)] -1 = F^iUn + \i - E(iu n )]. (7) 

k k 

Here F oc (z) is the scaled complimentary error function of a complex argument^ 

\ roo exp( — ?/ 2 ) 2 

Foo(z) = —= / dy = — isgn[lm(z)]\^Ke~ z erfc{— zsgn[Im(z)]z}. (8) 

y 7r J-oo z — y 

The dynamics of the (local) impurity problem is identical to the dynamics of the Ander- 
son impurity modelil'&B'ii and is determined by employing a weak-coupling conserving 
approximation for the local problem and satisfying the self-consistency relation in Eq. (0). 

It is important to note that since one does not a priori know the bare Green's func- 
tion Gq 1 in Eq. (||), one must iterate to determine a self-consistent solution for the Green's 
function of the infinite-dimensional lattice. This is done by performing self-consistent per- 
turbation theory for the self energy E[G] within a conserving approximation, and then 
determining the new local Green's function from the approximate self energy and Eq. ([7|). 
This process is iterated until convergence is achieved [the maximum variation of each G(iou n ) 
is less than one part in 10 8 which typically takes between 5 and 30 iterations]. 

Static two particle properties are also easily calculated since the irreducible vertex func- 
tion is locafi The static susceptibility for CDW order is given by 

X CW (q) = ^ E e^-^T fdr f dr'[(n ja (r)n k Ar')) - K.(r))(n^(r'))] 

ErEf w (q,^ m) ^) = TExr(q) , (9) 

mn mn 

at each ordering wavevector q. Dyson's equation for the two-particle Green's function 
becomesS0 

X C mn W (q) = Xl^rnn ~ ^ xKqjC^W > ( 10 ) 

P 

with the (local) irreducible vertex function in the CDW channel. 

The bare CDW susceptibility x°(q) i n Eq. ( |TUD is defined in terms of the single-particle 
Green's function 
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Xn(q) = -^E G n( k ) G n( k + q) 
iV k 



1 r°° . e-y 2 



dy- 



iuj n + fi-E n - X(q)y 



1-X 2 (q) 



yjl - X 2 (q) ^ n + /i-S n 

and all of the wavevector dependence is included in the scalaiHl X(q) = Ej=iCOsqj/d. 
The mapping q i— > X(q) is a many-to-one mapping that determines an equivalence class of 
wavevectors in the Brillouin zone. "General" wavevectors are all mapped to X = since 
cosqj can be thought of as a random number between —1 and 1 for "general" points in 
the Brillouin zone. Furthermore, all possible values of X (—1 < X < 1) can be labeled by 
a wavevector that lies on the diagonal of the first Brillouin zone extending from the zone 
center (X = 1) to the zone corner (X = — 1). The presence of incommensurate order in the 
attractive Holstein model is restricted to a very narrow region of parameter spaceii so only 
the "antiferromagnetic" point X = — 1 is considered for CDW order. The integral for x m {X) 
in Eq. fll!]) can then be performed analyticallyil = —1) = — G n /(iu n + fj, — E n ) . 

The irreducible vertex function is calculated perturbatively from the dressed Green's 

functions in a conserving approximation (see below). 

A similar procedure is used for the singlet s-wave SC channel. The corresponding defi- 
nitions are as follows: The static susceptibility in the superconducting channel is defined to 
be 

r/3 r/3 



x sc { 



iV -i-> t-> J J 



ET^f(q lWm)W „) = T^xSq) , (12) 

ran ran 

for superconducting pairs that carry momentum q; Dyson's equation becomes 

*22(q) = (q)<U - T]r^'(q)rSx P 5 „ c (q) > ( 13 ) 



with the corresponding irreducible vertex function for the SC channel; the bare pair-field 
susceptibility becomes 

xr(q)4E G »( k ) G -i(- k+( i) 

iv k 



-y 2 



dy~ ; ^ F t 



iuj- n -i + \i - S-„-i - X(q)y 



l-^ 2 (q) 



(14) 



yl - X 2 (q) iuJ n + ^ - S n - y 

with the special value Xn'(X = 1) = G n j{iio^ n ^\ + fi — S_ n _i) for the SC pair that carries 
no net momentum; and finally the irreducible vertex function is also determined in the 
conserving formalism (see below). 

At this point the transition temperature of the infinite-dimensional Holstein model is 
found by calculating the temperature at which the relevant susceptibility diverges (CDW or 
SC). 

Section II contains the comparison of QMC exact solutions to ME theory and the second- 
order conserving approximation for the Holstein model. Analytical expressions for the 
change in T c due to vertex corrections are given for the SC channel. Section III includes the 
application of conserving approximations to the attractive Hubbard model at half-filling. 
Truncated conserving approximations through fourth order are compared to the different 
fluctuation-exchange approximations and the exact QMC solutions. Conclusions are pre- 
sented in Section IV. 

II. HOLSTEIN MODEL 

There are two different types of approximations that are generally made for the electron- 
phonon interaction: the first method is a truncated conserving approximation that in- 
cludes all vertex corrections to a finite order and is valid for all values of the phonon 
frequencyElH'il'Ei; the second method is ME theory in which vertex corrections are ne- 
glected, but the phonon propagator is dressed to all ordersffiB These two methods are 
compared and contrasted in Figures 1 and 2. Figure 1 (a) shows the self energy for a con- 
serving approximation through second order. The self energy includes, respectively, the 
Hartree term (which is a constant and can be reabsorbed into the chemical potential), the 
Fock term, the second-order term that dresses the phonon propagator, and the lowest-order 
vertex correction. Figure 1 (b) displays the corresponding self-consistent equations for ME 
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theory: the self energy includes the Hartree term (which can once again be reabsorbed into 
the chemical potential) and the Fock term (which is calculated with the dressed phonon 
propag atorS). The dressed phonon propagator satisfies Dyson's equation [Figure 1 (c)]. 
To be more explicit, the self energy for the second-order conserving approximation is 

S n (cons.) = -g 2 Tj2 D n- r G r + g A T 2 ^[-2£> n _ r + D r _ s ]D n _ r G r G s G n - r+s , (15) 

r rs 

which includes the Fock diagram contribution and the two second-order contributions in 
Figure 1 (a). The bare phonon propagator D\ = Diiuji) in Eq. (|TJ|) is given by 
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A = - „ /na , , . , (16) 



for each (bosonic) Matsubara frequency u>i = 2lnT. On the other hand, in ME theory, the 
self energy satisfies 

W °^ J/(ff t ^) + l^ ; ^^T^G l+rGr (17) 
with TI; the phonon self energy (evaluated in the limit where vertex corrections are neglected). 

The self-consistency step involves determining a new local Green's function G n from the 
integral relation in Eq. (^) with the approximate self-energy of Eq. (15) or Eq. fllTf) . This 



process is repeated until the maximum deviation in the local Green's function is less than 
one part in 10 8 . 

Once the Green's functions and self energies have been determined, the irreducible vertex 
functions can be calculated for the CDW or SC channels. The vertices are simple in the ME 
theory: the irreducible vertex in the CDW channel satisfies [see Figure 2 (a)] 

T c m D n w (ME)=2U ; (18) 

and the irreducible vertex function in the SC channel satisfies [see Figure 2 (b)] 

r - (M£ ^- g2 M(^ + ^L) + n„,_„ ' < 19 > 

The phonon propagator in the CDW vertex is the bare propagator to avoid a double counting 
of diagrams. 
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The irreducible vertex functions acquire more structure in the second-order conserving 
approximation. In the CDW channel [see Figure 2 (c)] one must include both direct and 
exchange diagrams as well as the vertex corrections. The result is 

T c Ji w {cons.) =2U- 2U 2 TY,[G r G m . n+r + G r G m+n ^][ ^ 2 ] 2 

Si 2 il 2 Q 2 Si 2 
tt _ oTT 2 r\^rr f 1 

W + ^-n V UrUm ~ n+r W + U 2 m _ n [ W + U 2 m _ n W + Ulj 

il 2 Vt 2 

+ U 2 T £ G r G m+n _ r Q2 - ^- fi2 - ^- ■ (20) 

Note that the vertex corrections (arising from the first-order exchange diagrams) modify the 
interaction in the CDW channel so that it properly interpolates between the zero frequency 
limit T CDW — > 2U and the infinite frequency limit r — > [/. At an intermediate frequency, 
the CDW interaction strength has a complicated temperature dependence. In the SC channel 
[see Figure 2 (d)] one finds 

Q 2 

T s m c n (cons.) = U 2 

Si T" Ul m _ n 

2 il 2 il 2 il 2 il 2 

+ [/ T G r G m ~ n+r — 2^ — [2 



il 2 + Wm_ n i! 2 + Wm_n il 2 + UJ 2 ~r ^ 2 + w i+r+l 

si 2 Q 2 

- U 2 T ^ G r G m+n+r+1 . (21) 

As the transition temperature (to a CDW-ordered state or a SC-ordered state) is ap- 
proached from above, the susceptibility (in the relevant channel) diverges. Therefore, one 
can determine the transition temperature by finding the temperature where the scattering 
matrix (in the relevant channel) 

T m n XT 'mnXn ' (22) 

has unit eigenvalue^. In general, the eigenvector corresponding to the maximum eigenvalue 
of the scattering matrix is symmetric with respect to a change in sign of the Matsubara 
frequency. 

At half filling the Holstein model interaction is particle-hole symmetric, so the Green's 
functions and self energies are purely imaginary and the vertices are real. The self energy 
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can be expressed by T,(iu n ) = iou n Z(iuj n ), with Z(iui n ) the renormalization function for the 
self energy. At half filling, both = ~ 1) anc ^ Xm(X = 1) are also even functions of 

the Matsubara frequency, so the only contribution of the irreducible vertex function to the 
eigenvalue of the scattering matrix comes from the even Matsubara frequency component 

In order to judge the accuracy of these approximate methods for the electronic self energy 
and the irreducible vertex functions, it is necessary to compare them to the exact results. 
The best way to do this would be to directly compare the perturbative results to the exact 
QMC results. Unfortunately, there is no available QMC data to do this. However, it has 
been well established, that the iterated perturbation theory (IPT) of Georges and KotliaiS 
yields essentially exact results for the electronic self energy of the Hubbard model (by direct 
comparison with the QMC results@) as long as the system is at half filling. The IPT is 
identical to the second-order conserving approximation, except that the perturbation theory 
is strictly truncated to second order in U . 

A comparison of the approximations to the IPT results for the electronic self energy 
and a comparison of the approximations for one column of the irreducible vertex function 
in the CDW channel is made in Figures 3 and 4 for two different interaction strengths at 
half filling. The phonon frequency is set to be approximately one-eighth of the effective 
bandwidth (Q/t* = 0.5) as was done in the QMC solutions!. The energy cutoff is set to 
include 256 positive Matsubara frequencies for the perturbative approximations. 

At weak coupling (g = 0.4, Figure 3), the second-order conserving approximation clearly 
provides a more accurate approximation to the electron self energy (under the assumption 
that the IPT is essentially exact). One expects the exact irreducible vertex function to be a 
frequency dependent interaction, so the second-order conserving approximation is probably 
more accurate here too [the CDW vertex for the ME theory has no frequency dependence 
as shown in Eq. flIS|)1. 

An underestimation of the self energy causes an overestimation of the transition tempera- 
ture and vice versa. Similarly, an underestimation of the magnitude of the irreducible vertex 
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will cause an underestimation of T c (and vice versa). Since the ME theory overestimates 
the self energy and overestimates the CDW vertex (because there is no weakening of the 
CDW vertex at small frequency transfers) , it is difficult to predict whether or not ME theory 
will overestimate T c . In the same fashion we do not know whether or not the second-order 
conserving approximation will overestimate T c since the magnitude of the exact vertex is 
not known. 

As the coupling strength is increased to the point where a double-well structure began 
to develop in the effective phonon potential of the QMC simulations^ (g = 0.5, Figure 4) 
one can see strong-coupling effects begin to enter. Surprisingly, ME theory is becoming a 
more accurate approximation to the self energy here. Stated in other words, as the coupling 
strength increases, the effect of vertex corrections is reducedS. Unfortunately, the self- 
consistent equations for the ME theory become unstable to an iterative solution as the 
coupling strength is increased further. 

At half filling, the Holstein model always has a transition to a CDW-ordered phase at 
q = (jr, 7T, 7T...) (X = —1). The transition temperature to this commensurate CDW is 
plotted in Figure 5 as a function of the interaction strength. The second-order conserv- 
ing approximation is compared to ME theory and the QMC simulationslil. The conserving 
approximation is much more accurate at weak coupling@0Jii (ME theory predicts a tran- 
sition temperature that is an order of magnitude higher than the QMC and conserving 
approximation results at the lowest value of the coupling considered) because the inclusion 
of the first-order exchange diagrams produces the correct interaction and the inclusion of the 
second-order terms produces the correct prefactor. However, ME theory does display the 
proper qualitative behavior of developing a peak in T c as the interaction strength increases. 
This feature is not reproduced by the truncated conserving approximation. 

As the system is doped away from half-filling, the CDW instability remains locked at the 
commensurate point (X = —1) until it gives way to a SC instability (incommensurate order 
may appear in a very narrow region of phase space near the CDW-SC phase boundary0il 
but is neglected here). In Figure 6, the phase diagram of the Holstein model is plotted for 
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two values of the interaction strength (g = 0.4, g = 0.5). The weak-coupling QMC data 
(g = 0.4) are reproduced most accurately by the second-order conserving approximation, 
as expected from the comparison of the self energy and the vertices in Figure 3. The SC 
transition is reproduced remarkably well, because the underestimation of the self energy 
[Figure 3 (a)] must be compensated by an underestimation of the SC vertex. The critical 
concentration for the CDW-SC phase boundary is also more accurately determined by the 
conserving approximation. Note that the difference between the SC transition temperature 
calculated with ME theory and the second-order conserving approximation explicitly shows 
the lowest-order effect of vertex corrections. The vertex corrections lower T c by about a 
factor of two at the phase boundary, but are reduced as the doping increases. 

At the stronger coupling strength (g = 0.5) ME theory reproduces the CDW transition 
temperature very accurately at half filling, but not the doping dependence of T c . It does 
manage to reproduce the transition temperature in the SC sector quite well, but the second- 
order conserving approximation is superior at determining the CDW-SC phase boundary. 
Clearly both approximation methods are inadequate at this large a value of the coupling 
strength. 

Up to this point we have concentrated on one value of the phonon frequency and have 
compared the numerical solution of the self-consistent perturbation theory with the nu- 
merically exact QMC solutions. In the limit of weak coupling (\U\ — ► 0), the transition 
temperature approaches zero (T c — > 0) and the leading behavior of T c can be determined 




ion. 

We will concentrate on the SC channel only, b 



ecause the analytical 



techniques are not as accurate for the CDW channel. 



In the limit T — > 0, the self energy satisfies 



lim [iu n + n — T:(iu n )] = \i + iui n Z 



(23) 



where the renormalized chemical potential is 



fi = fx — lim KeL(iuj n ) = \i — Up e + 0(U 2 ) 



(24) 
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and the renormalization function is 

Z = Z(0) = 1 + \U\ f iy »^ + "> Y'"" + " } (n " yf + ■ < 25 > 

Here p(y) = exp(— y 2 )/ yjir is the noninteracting density of states (DOS) in infinite dimen- 
sions. The irreducible vertex in the SC channel [Eq. fl2"ID l becomes 

r s r f n = - TW -f^\U\il + 2\U\h]+U 2 I 2 + 0(U 3 ) , (26) 

1 1 + UJ m _ n 

with I\ and I 2 two smooth temperature-dependent integrals that can be approximated by 
their zero-temperature limit: 

1 r°° f] 4 

1 ^-tJ-J^ + ^w^? + ° {u) ■ (28) 

The bare susceptibility becomes 

x « (X = l) = - lmF ^f + » + 0(U>) , (29) 



for the electron pairs that carry no net momentum. 

In the square-well approximation^'!, the smooth temperature dependence of the SC 
vertex is replaced by a sharp cutoff at a characteristic frequency u c : 

Q 2 

Q2 2 > 0(lj c - \u m \)0(u c - \uj n \) , (30) 

1 L -+- Ul m _ n 

with 6(x) the unit step function. The scattering matrix [in Eq. fl22|) 1 then becomes 



T mn = -T[6{u c - \u m \)6{u c - K|)|C/|(1 + 21171/0 - ^/ 2 ] ImF °°^ + ^ . (31) 



The SC transition temperature is now determined by solving the matrix eigenvalue equation 

The eigenvector <f) m , can be chosen to be of the form (fi m = 1 + a9(u c — \ui m \) in the 
square-well approximation, so that the matrix eigenvalue equation is reduced to two coupled 
algebraic equations 
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a = pitM {1 + mh){1 + a)R ? 

l = _p_itEl h{aR + s) f (32) 

with R and £ denned by 

5- = T InxF^ioyg + g) . . 

p(/i) n ±^ w„ 

The infinite summation over Matsubara frequencies can be performed in the standard 
fashion! to yield 

S=lr ^ + /i) tanh. V 



2 J-oo y p(n) 2ZT C 



— + / — (tanhy-l + , 35 

ZT C Jo y 2p(p) 

while the truncated summation can be expressed in an integral form0 



R = S-^-[ —ImF^iiyZ + n) , 

Itx Juc y 

= s _i_ r° dy P (y + g) + p(-y + g) tan _ 1 

7r jo ?/ 2p(p) zu c 

if the transition temperature is much less than the cutoff frequency (T c << uj c ). 
The coupled algebraic equations [in Eq. (|32|)1 are then solved by 

-^_ = [l + |r/|(2/ a -/ 2 )]5+ (l + 2|l/|ii)(i2 -S) , (37) 

to order \U\. The limiting form of the transition temperature is now found by substituting 
Eqs. (|35| ) and for i? and 5 into Eq. ( p7|) and solving for T c . The result is 

T c GXp( - ) fphononf electron f DOS f vertex j (38) 

which includes the interaction term and the constant prefactors. The constant terms arise 
from the phonon self energy, the electron self energy, the nonconstant DOS, and the vertex 
corrections. 
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The phonon self energy correction is 

f phonon = exp[-— -^-yRe J dyF^(iy + /i)] = exp[V2 J — erf(w + v2//)] , (39) 

and is independent of the phonon frequency. This correction factor is normally included in 
the definition of the electron-phonon interaction strength A 

l/T I poo 

\ = p ^)\U\[l- l -^Re dyFl(iy + fJL)] , (40) 

71 J-oo 

of the ME theory formalism. 

The electron self energy factor satisfies 

f°° 2 Q 

/electron = 6Xp[- J ^ dyt~ V COSh(2/Jg/) + ] , (41) 

which approaches the standard ME theory result of j "electron - > e -1 as Q — > and approaches 
the standard Hubbard model result of f electron — > 1 as Q — > oo. 
The nonconstant DOS factor is 

1 f°° dv 2 2 v 

f D os = 7T ex P / — t tanh y - 1 + e ~ V cosh(2/iy) (1 - - tan" 1 JL)] (42) 

2 JO y 7T LO c 

which depends on the square-well cutoff frequency u c . In the limit u c — > 0, /dos approaches 
the ME theory result of 1.14o; c , whereas in the limit u; c — > oo, foos approaches the Hubbard 



model resultS 

lim f DOS = 0.85 exp[2 [°° ^-e~ y2 sinh 2 \iy\ = 0.85 exp[v^F f dye*ed(y)\ . (43) 

Jo y Jo 



Note that once the phonon energy scale is larger than the electronic energy scale, it is the 
bandstructure (not u c ) that determines the DOS prefactor. 
Finally, the vertex correction factor becomes 

1 roc Q 2 I tt 2 

/^ = exp[^Rey_^(z, + ^)^ T ^(l + -^ T -^)] , (44) 

which approaches 1 as Q — > 0. In the high-frequency limit, the vertex corrections cancel the 
phonon self energy corrections and yield 
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lim fphononfvertex = exp[ — \= ( — ^j^-erf(w + V2fj,)} , (45) 

which reproduces van Dongen's resultEl fphononfvertex = exp[— \/2hi(l + \/2)] at half filling 
(/i = 0). In the infinite frequency limit (attractive Hubbard model), the vertex corrections 
always reduce T c . 

The effect of vertex corrections upon the superconducting transition temperature in the 
weak-coupling limit are displayed in Figure 7 (a). The vertex correction factor, f ver tex = 
T c (vertex) /T c (no vertex), is plotted against the phonon frequency for eight different electron 
concentrations. At half filling, the vertex corrections sharply reduce T c , so that T c calculated 
with vertex corrections is a factor of two lower than T c calculated without vertex corrections 
at Q = 0.13t* (or, since the effective electronic bandwidth is approximately W = 4t*, when 
Q/W = 0.03). Therefore, vertex corrections should play an important role in Bai-^K^BiOs 
where tt/W = 0.02. 

As the system is doped away from half filling the effect of the vertex corrections is reduced 
(as was already seen in Figure 6), until a critical electron concentration (p c « 0.2) is reached 
where the vertex corrections initially cause an enhancement to T c . This enhancement occurs 
because the electronic Green's functions have a larger real part than imaginary part which 
causes the integrand in Eq. flli]) to change sign for small y. This enhancement will not 
be seen in standard ME theory with a constant DOS, because the Green's functions are 
chosen to be purely imaginary in that case. As the phonon frequency is increased to a large 
enough value, the vertex corrections will once again reduce T c , because they always cause a 
reduction in the limit Q — > oo [see Eq. (§^)]. 

The square-well cutoff frequency, u> c , should vanish as the phonon frequency vanishes, 
and should become infinite as the phonon frequency becomes infinite. The cutoff frequency 
is chosen to be three-fifths of the phonon frequency (ou c = 0.6f2), so that the proper limiting 
behavior is attainediiS as fl — ► and £1 — > oo. The prefactor to the SC transition 
temperature in Eq. ( RSI ) is plotted in Figure 7 (b) for eight different electron concentrations. 
Note that there is an optimal phonon frequency where the SC response is maximal, that 
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shifts from a low frequency at half filling to higher frequencies as the electron concentration 
is reduced. The optimal phonon frequency is usually smaller than the effective bandwidth 
(W = At*). 

In general, one must expand the free energy thru second orderS^ to properly determine 
T c in the limit \U\ — ^ 0. The miracle of ME theory is that a first-order calculation with 
dressed phonons properly determines T c in the limit Q — > 0. In the CDW channel, the vertex 
corrections modify both the interaction strength and the prefactor, while in the SC channel, 
only the prefactor is modified. It is this robustness of the SC channel to the effects of vertex 
corrections that explains the success of ME theory for low temperature superconductors. 



III. HUBBARD MODEL 

The Hubbard model in Eq. (0) is the infinite-frequency limit (Q — > oo) of the Holstein 
model. The Hubbard model has an electron-electron interaction that only occurs between 
electrons with opposite spins. This happens because of the cancelation of the direct and 
exchange diagrams which causes all electron-electron interactions between like-spin parti- 
cles to vanish. The perturbation theory becomes much simpler in the Hubbard model case, 
because of this reduction of diagrams, and can be performed to higher order. Here the trun- 
cated conserving approximation will be carried out to fourth order, and will be compared 
to the fluctuation-exchange (FLEX) approximationll! to determine the best way to approxi- 
mate the Hubbard model in the infinite-dimensional limit. Previous work has concentrated 
on second-order conserving approximations^'^, third-order conserving approximations^, or 
the FLEX approximation!!. 

One expects that a truncated approximation will be superior to an infinite summation of 
random-phase approximation (RPA) bubbles and particle-hole and particle-particle ladders 
because the many-body problem reduces to a self-consistently embedded Anderson impurity 
model, and the analysis of Yamada@ has shown that the total fourth-order corrections to 
the self energy are an order of magnitude smaller and opposite in sign to the fourth-order 
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contribution of the FLEX approximation. The irreducible vertex functions should have 
similar effects, but have not yet been analyzed in detail. 

The diagrammatic expansion for self energy (in a conserving approximation) of the Hub- 
bard model is given in Figure 8. The first line includes the first-order Hartree contribu- 
tion (which is a constant that can be absorbed into a renormalized chemical potential), 
the second-order contribution and the third-order particle-hole and particle-particle lad- 
ders. The second line contains the fourth-order contributions from the RPA bubbles and 
the particle-hole and particle-particle ladders. The third and fourth lines include all of 
the remaining fourth-order diagrams (the inclusion of the second-order self energy into the 
dressed Green's function of the second-order diagram produces another fourth-order con- 
tribution to the self energy, but this is automatically included in the self-consistency step 
of the conserving approximation). The FLEX approximation consists of the summation of 
all RPA bubbles, particle-hole ladders, and particle-particle ladders. The self-energy has 
already been determined on the real axis by Menge and Muller-Hartmannil. The FLEX 
approximation for the self energy includes all contributions thru third order in U (the first 
line of Figure 8) but only a partial contribution of the fourth-order and higher-order terms 
(the second line of Figure 8 plus the higher-order terms). An explicit formula for the elec- 
tronic self energy of the Hubbard model thru fourth order is given in the appendix. The 
corresponding formula for the FLEX approximation has been given befordll'Il. 

The irreducible vertex functions are too cumbersome to represent diagrammatically, but 
an explicit formula for the CDW vertex is given in the appendix. The corresponding for- 
mula for the FLEX approximation has already been giverJll. An additional simplification 
is usually made for the FLEX approximation that neglects a large class of diagrams for 
the irreducible vertices (the so-called Aslamazov-Larkin diagram thereby including only 
those contributions to the irreducible vertex function that can be represented by functions 
of the bare particle-hole or bare particle-particle susceptibilities^!. This simplified FLEX 
approximation will be denoted FLEX*. 

Since the Hubbard-model interaction is particle-hole symmetric, the half-filled band cor- 
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responds to fi = 0, and the Green's functions are purely imaginary. The odd-order contribu- 
tions to the self energy all vanish and each of the fourth-order contributions on a given line 
in Figure 8 are identical^ (see the appendix). Furthermore, it can be shown that the most 
complicated contributions to the irreducible vertex function in the CDW channel are odd 
under a change in the sign of the Matsubara frequency, and can be neglected in calculating 
the maximum eigenvalue of the scattering matrix, because only the even component of the 
irreducible vertex function enters (see the appendix for details). 

Note that since the self-energy is an odd function of U at half-filling, but the irreducible 
vertex function contains both even and odd powers of U, the only difference between a 
truncated conserving approximation of order 2n and of order 2n + 1 is that the irreducible 
vertex function is larger for the odd-order approximation. Therefore, we expect that an even 
order approximation will underestimate the transition temperature (in weak-coupling) and 
an odd-order approximation will overestimate T c . 

A comparison of the different approximation schemes is given in Figures 9 and 10 for two 
different values of U. The second-order, third-order and FLEX approximations all employ 
an energy cutoff of 256 positive Matsubara frequencies; the fourth-order approximation uses 
64 positive Matsubara frequencies. In Figure 9 (a) the self-energy renormalization function 
is plotted for the three different approximations at U = —t* and compared to the essen- 
tially exact IPTS. Note that the fourth-order approximation virtually reproduces the IPT 
results, but that the FLEX approximation grossly overestimates the self-energy even though 
the coupling strength is not too large. In Figure 9 (b) the even component of one row 
of the irreducible vertex function for the CDW channel is compared for U = —t*. All of 
the truncated conserving approximations are in reasonable agreement with each other; the 
FLEX approximation has the smallest magnitude at low Matsubara frequency. The simpli- 
fied FLEX* grossly overestimates the magnitude of the vertex (in fact the FLEX* approx- 
imation produces the wrong qualitative behavior of the vertex). In general, the transition 
temperature calculated with the simplified FLEX* will be a more accurate approximation 
to the exact T c than that calculated with the full FLEX, because the overestimation of the 
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self-energy will be compensated for by the overestimation of the vertex in FLEX*. These 
results are similar to what found for the repulsive Anderson impurity model. 

As the coupling strength is increased to U — —2t*, the FLEX approximation becomes a 
more accurate approximation for the self energy than the truncated conserving approxima- 
tions [see Figure 10 (a)]. Clearly, the truncated conserving approximation must be performed 
to a high order to accurately reproduce the self energy in the strongly correlated regime. The 
irreducible vertex function in the CDW channel is plotted in Figure 10 (b). The different 
approximations no longer agree well with each other indicating that the perturbation theory 
is breaking down. 

The transition temperature for the CDW transition at half-filling in the attractive Hub- 
bard model is plotted in Figure 11 as a function of U. The truncated conserving approx- 
imations are doing quite well in the weak-coupling regime. The odd-order approximation 
overestimates T c , while the even-order approximations underestimate T c . The even-order 
approximations tend to be more accurate over a wider range of U than the odd-order ap- 
proximation, but neither approximation properly reproduces the turnover in T c as a function 
of U as seen in the QMC simulationsil. Note also that all truncated approximations agree 
in the limit U — > 0, but that a first-order (RPA) calculation will be off by a factor of three 
in the weak-coupling limitfflS. The FLEX approximation does have the correct quali- 
tative behavior of developing a peak in T c as a function of U but the peak position and 
peak height are off by about an order of magnitude. The simplified FLEX* is, in general, 
a more accurate approximation than the full FLEX, but becomes unstable if U is increased 
too far. The FLEX* does not agree as well with the QMC calculations (or with the other 
approximations) in the weak-coupling limit because the irreducible vertex function does not 
include all of the third-order contributions. 

The reason why the truncated (even-order) conserving approximations do not approx- 
imate the self energy (or the vertex) too accurately at moderate coupling, but are good 
approximations for the transition temperature at moderate coupling is most likely due to 
a cancellation of the effect of an underestimation of the self energy by an underestima- 
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tion of the vertex in the calculation of T c . This probably explains why the T c curves do 
not turnover for the truncated approximations as well: the turnover must be arising from 
self-energy effects that are being underestimated here. 

In summary, the truncated approximations tend to give better numerical agreement than 
an approximation that tries to sum an infinite series of diagrams (such as the FLEX). One 
must go to very high order to see a peak develop in T c as a function of U and to have good 
quantitative agreement with the QMC results. It will be interesting to see if the removal of 
the infinite summation of diagrams in the conserving approximation produces an even better 
agreement with the QMC results (as was found for the Anderson impurity modefl). The 
approximation will no longer be a conserving one, and will need to be generalized to move off 
of half-filling, but should be even more accurate. The fluctuation-exchange approximation 
seems to be a poor approximation, and should not be tried for the Holstein model, rather 
one should concentrate on generalizing Yamada's analysis for instantaneous interactions to 
one for retarded interactions to see whether or not one can improve upon the accuracy in 
that case too. 



IV. CONCLUSIONS 

Vertex corrections can be systematically incorporated into a weak-coupling theory of 
electron-phonon interactions. Expansions must be performed to second order in the effective 
electron-electron interaction in order to produce the correct behavior in the weak-coupling 
limitffl. The miracle of Migdal-Eliashberg theoryiS is that a first-order calculation suffices 
(with dressed phonon propagators) in the small-phonon- frequency limit. Vertex corrections 
enter to lowest order in the CDW channel, modifying the interaction strength. They enter 
to higher order in the SC channel, and merely modify the prefactor of the weak-coupling T c 
equation. This robustness of the SC channel to the effects of vertex corrections explains its 
remarkable success for low temperature superconductors. Nevertheless, the effect of vertex 
corrections should be strong enough to be observable in materials such as Ba 1 _ :r K :r Bi03 and 
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the doped fullerenes. It is possible that the effects of vertex corrections can even be detected 
in certain low-temperature superconductors such as Pb. 

In general, vertex corrections will reduce the transition temperature, however there is a 
small parameter regime at low electron density where the vertex corrections actually cause 
an enhancement to the superconducting transition temperature. This occurs because the 
real parts of the Green's functions are larger than the imaginary parts for small imaginary 
frequency and low electron concentration. At high enough phonon frequency, or large enough 
electron density, the vertex corrections will lower T c . 

Truncated conserving approximations appear to be better approximations than infinite 
summation schemes such as the fluctuation-exchange approximation^. The electronic self 
energy, the irreducible vertex functions, and the transition temperatures all appear to be 
better approximated by a truncated conserving approximation. The qualitative feature of the 
development of a peak in the transition temperature as a function of the interaction strength 
is, however, not reproduced by a truncated conserving approximation. Perhaps a completely 
truncated approximation (that is no longer conserving) will do even better at approximating 
properties of interacting electronic systems in infinite dimensions. Yamada^0 found this 
to be so for the Anderson impurity model, and his techniques have been applied in infinite 
dimensionsBil to second order in U. What is needed is a way to generalize Yamada's work 
off of half filling for both the self energy and the irreducible vertex functions. Work in this 
direction is currently in progress. 

In conclusion, a weak-coupling conserving approximation has been carried out for the 
attractive Holstein and Hubbard models that includes all effects of vertex corrections and 
nonconstant density of states. Agreement with the exact solutions is found to be excellent 
at weak-coupling, but the qualitative feature of developing a peak in T c as a function of the 
interaction strength is not reproduced. From this standpoint, a weak-coupling theory is much 
more difficult to control than a strong-coupling theory (perturbation theory in the kinetic 
energy). Analytic expressions for T c in the SC channel have been explicitly derived, and 
they indicate that vertex corrections may be observable for some classes of low temperature 
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superconductors. 

Future work will include an examination of the ordered phase, a study of the effects 
of Coulomb repulsion, and a real materials calculation to look for the effects of vertex 
corrections in low-temperature superconductors. 
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APPENDIX A: APPENDIX 

The explicit formulas for the electronic self energy thru fourth order and for the irre- 
ducible vertex function in the CDW channel are given here. The diagrammatic expansion 
for the self energy is shown in Figure 8. The first two lines include the contributions from 
the FLEX approximation (truncated at fourth order) while the last two lines include the 
extra fourth-order contributions. The self energy is then expanded as 

E(iu n ) = S n = EX (4) + EU4) , (AI) 

with E ELEX (4) the contributions to the self energy included in the FLEX approximation^! 
(but truncated to fourth order) and E^(4) the additional contributions to fourth order. The 
truncated FLEX contributions ar 

e flex (4) =TU <£ Gn _ t [Ux f + u^h2 + 2U *tf*\ +TUJ2 GL*-i+, [U 2 X f 2 - U 3 x7 3 ] , 
i i 

(A2) 
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thru fourth order. Here the bare particle-hole and particle-particle susceptibilities are 

X ph {i^i) = xf = -T^G r G r+ i , X w N = Xr = rE^ G -i + i • (A3) 

r r 

The additional fourth-order contributions to the self energy are 

^ ) = T 2 U 4 ^2 Gn-iGn-i+iiGn+vx^X^* 1 + 2T 2 U 4 ^ G ^ n ^ 1+i G ^ n - 1+ i-i'G n+ i>x PP xln+i-i+i' 



T 3 U 4 G n -iG r G r+ iG r+ i + iiG r+ i'X P ^ — T 3 U 4 ^ G n ^iG n -iiG r G r -i + iiG r+ i'X P 



rlV rlV 



_l_ ^3^4 ^2 G n ^iG n -i'G r G r+ iG r+ iix P n- m -i-ii 

rlV 

At half filling the Green's functions are purely imaginary and satisfy GL n -i = —G n . 
Therefore, the particle-hole and particle-particle susceptibilities in Eq. (|A3|) are equal. It 
is easy to show that the two third-order contributions to S FLEX (4) vanish in this case, and 
that the three fourth-order contributions are equal. Similarly, the three terms that involve 
a double summation in Eq. are equal, and so are the three triple summation terms as 
shown by Yamadall. 

The irreducible vertex function in the CDW channel can also be determined. The vertex 
function is broken up into its FLEX contributions and its additional fourth order contribu- 
tions 

r^ w = r^ EX (4) + C„(4) . (as) 

The FLEX contributions thru fourth order areEl 

r!nt EX ( 4 ) = Z7 + ^Xm-n[2 + UXm-n + ^^Xm-n] ~ ^Xm+n+l^ — ^Xm+n+l + ^Xm+n+l] 

— 2TU 2 ^ G r G n - m+r [U Xm-r + 3U 2 Xm-r] ~~ 2TJJ 2 ^ G r G m+n - r [U Xm-r + 3U 2 Xm-r] 

r r 

— 2TU 2 G r G m - n+r [—2U Xm+r+l + 3t/ 2 Xm+ r +i] , 

r 

and the additional fourth-order contributions are 

-p' ( a \ — ottt a X^ r r \v ph v pp -i_ v p/l v ph -i- v ph v pp 1 

mn\ ) u / , KJr r^ J n~rn+r[A.m—rA.n+r+l ' Am—nXm—r T Am-nAn+r+lJ 

r 

, rpjjA sr^ r r \O^P h ^PP I v?' 1 ^P h 1 

~r -L <-J / t ^rum+n-r [■'Ara-rAm+n+l ' Am-rAn-rJ 
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_i_ /It 2 77 4 r r r r w ph -u ^' ph i 

' ^ tJ - u / j IJ r' J s u n-r+s u n-m+s [Am-r ~T~ Xn—ri 
rs 

— 9T 2 77 4 V"* r 1 r* r 1 rr* -v ph 4- c v ph 1 

z, - t u / j ^r^ s yjc m+n—s\ KJ n+i — sXm—r ~r v - r m+r—sXn—r] 
rs 

_ n T 2 T ji n n n n \ v ph _i_ v P h l 

z, - t ^ / y u r u s' J -m+r+s lJ ii-m+s[Am-r ~~r Xn—ri 
rs 

T U ^ ' GfGg [(jn-m^-rGn-m^s G m -\- n — r G m -\- n — s ] Xr— s 
rs 

2X 1 tT" ^ ] G r G s G m + r — s G n + r — s \"X.m+r+l Xn+r+l] T U ^ ] G r G S G n — m + r G m — n+sX^r+s 



rs rs 

2T U ^ ' G r G s \G m J r ^ — s G n — r+sXm+f+l ^n+i — s^m-r+sXn+r+l] 
rs 



rst 

_l_ ^3^-4 G r G s G t G r+s - t G 

n+s—t\G m —r+t ^m+r-f] 



rst 



Note that the simplified FLEX* approximation does not include any of the terms that involve 
explicit summations over Matsubara frequency [the second and third lines in Eq. (|A6|) ]. 

By making the transformations t — > r + s — t, r s, and s — > r in the triple summation 
terms in Eq. (|A7|) and using the symmetry at half filling G- n -i = — G n , one can demonstrate 
that the triple summation terms are odd under n — > — n — 1, and do not contribute to the 
eigenvalue of the scattering matrix if the eigenvector is even under n — > — n — 1. Therefore, 
the triple-summation terms in Eq. ( |A7| ) may be neglected (this result has been explicitly 
tested by calculating the eigenvalue of the scattering matrix with and without the triple- 
summation terms and there was no effect on the eigenvalue at half filling). 
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FIGURES 

FIG. 1. Dyson equations for the self energy of the Holstein model. The thick solid lines denote 
dressed (electronic) Green's function and the thin wavy lines denote the phonon propagator. The 
self energy (expanded out to second order in a conserving approximation) is depicted in (a) and 
includes the Hartree and Fock contributions, the second-order dressing of the phonon line and the 
lowest-order vertex correction. The self-consistent equation for the self energy in Migdal-Eliashberg 
theory is shown in (b). The thick wavy line is the dressed phonon propagator which satisfies the 
Dyson equation in (c). 

FIG. 2. The irreducible vertex functions in the CDW and SC channels. The CDW irreducible 
vertex function for ME theory is shown in (a). Note that the phonon propagator is bare in (a). The 
SC irreducible vertex function for ME theory appears in (b). Note that the phonon propagator is 
dressed here. The CDW irreducible vertex function for the second-order conserving approximation 
is shown in (c). Note that the vertex corrections (exchange diagrams) modify the interaction 
to lowest order in the CDW channel. The SC irreducible vertex function for the second-order 
conserving approximation is shown in (d) . The vertex corrections first enter at second order in the 
SC channel. 

FIG. 3. Comparison of the ME theory (solid line) to the second-order conserving approxima- 
tion (dashed line) for the Holstein model at half filling with phonon frequency 0, = 0.5t*, interaction 
strength g = OAt*, and temperature T = t*/lG. This example is generic for the weak-coupling 
limit. In (a) the self energy renormalization function Z(iu n ) — 1 is plotted against the Matsubara 
frequency and compared to the IPT (solid dots). In (b) the symmetric combination of the first 
column of the irreducible vertex function in the CDW channel is shown. Note that the second-order 
conserving approximation is clearly superior to ME theory in the limit of weak coupling. 
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FIG. 4. Comparison of the ME theory (solid line) to the second-order conserving approximation 
(dashed line) for the Holstein model at half filling with phonon frequency Q = OM*, interaction 
strength g = 0.5t*, and temperature T = t*/9. This example is generic for the transition region to 
the strong-coupling limit. The self-energy renormalization function (a) and the irreducible vertex 
function in the CDW channel (b) are both pictured. In (a) the self energy is compared to the IPT 
(solid dots). Note that in the limit where the strong-coupling effects begin to manifest themselves, 
the ME theory is becoming a more accurate approximation, or, put in other words, the total effect 
of vertex corrections is reduced as the interaction strength increases. 

FIG. 5. Transition temperature to the CDW-ordered state at half filling in the Holstein model 
at an intermediate phonon frequency (f2 = 0.5t*). The ME theory (solid line) is compared to the 
second-order conserving approximation (dashed line) and the QMC results (solid dots). Note that 
the vertex corrections are very important in the CDW channel and that only the second-order 
conserving approximation produces the correct result in the weak-coupling limit. ME theory does, 
however, display the correct qualitative behavior of developing a peak in T c as a function of inter- 
action strength. 

FIG. 6. Phase diagram of the Holstein model with Q = 0.5t* at two different coupling strengths 
(g = 0.4, 0.5). The solid dots are the QMC solutions with CDW order, and the open triangles are 
the QMC results with SC order. The kinks in the solid (ME) and dashed (second-order conserving 
approximation) lines occur at the CDW-SC phase boundaries. In the weak-coupling limit (g = 0.4) 
the second-order conserving approximation is superior to the ME theory and is quite accurate 
for the SC transition. The difference between the ME results and the second-order conserving 
approximation show explicitly the lowest-order effects of vertex corrections upon the SC transition 
temperature. The effect of vertex corrections is reduced as the filling is reduced. ME theory is 
quantitatively more accurate in determining T c for the stronger coupling strength (g = 0.5) but the 
second-order conserving approximation is superior in determining the CDW-SC phase boundary. 
Clearly both approximations are failing at such a large value of the interaction strength. 
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FIG. 7. Plot of the lowest-order effect of vertex corrections on the SC transition temperature in 
the weak-coupling limit. In (a) the renormalization factor for the transition temperature calculated 
with vertex corrections divided by the transition temperature calculated without vertex corrections 
is plotted as a function of the phonon frequency. Eight different values of the electron density are 
plotted (p e = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,1.0). Note that the vertex corrections always reduce 
T c in the high phonon frequency limit, but that vertex corrections initially enhance T c at low 
electron density (p e < 0.2). At half-filling (where the Green's functions are purely imaginary) 
vertex corrections will reduce T c by a factor of two when Cl/W = 0.03. In (b) the prefactor of 
the weak-coupling T c formula is plotted against the phonon frequency for the same eight values of 
the electron density. The cutoff frequency has been chosen to satisfy u) c = 0.6S7. Note that there 
is always an optimal phonon frequency where the SC response will be the largest, and that this 
optimal phonon frequency increases as the electron concentration decreases. 

FIG. 8. Self energy diagrams for the Hubbard model thru fourth order. The thick solid lines 
denote the dressed electronic Green's functions, and the thin dotted lines are the Coulomb inter- 
action. The first two lines contain all of the FLEX contributions truncated to fourth order. The 
last two lines are the remaining fourth-order diagrams. At half filling the odd-order contributions 
to £ vanish, and each of the three fourth-order contributions on the same line yield the same 
contribution to E. 

FIG. 9. Comparison of the different conserving approximations for the Hubbard model at 
half filling in the limit of weak-coupling (U = —t*,T = t*/20). The second-order (dashed line), 
third-order (dotted line), and fourth-order (solid line) conserving approximations are compared to 
the full FLEX (chain dotted line) and the simplified FLEX* (chain triple dotted line). In (a) the 
self-energy renormalization function is plotted against Matsubara frequency and compared to the 
IPT (solid dots). In (b) the even component of the first column of the irreducible vertex function in 
the CDW channel is plotted. Clearly the fourth-order approximation is the best approximation in 
this limit. The FLEX approximation grossly overestimates the self energy. The simplfied FLEX* 
compensates for this by overestimating the vertex to produce a more accurate value for T c . 
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FIG. 10. Same as in Figure 9, but with a stronger value of the coupling (U = —2t*,T = t*/8). 
In this limit the FLEX approximation is superior for the self energy, but the vertex does not appear 
to be reproduced accurately by any approximation. 

FIG. 11. Transition temperature to the CDW-ordered state in the Hubbard model at half filling. 
The second-order (solid line), third-order (dotted line), and fourth-order (solid line) conserving 
approximations are compared to the full FLEX (chain dotted line), the simplified FLEX* (chain 
triple dotted line), and the QMC results (solid dots). Note that the odd order approximations 
overestimate T c , the even-order approximations underestimate T c , and that one has to go to very 
high order to reproduce the peak in the transition temperature as a function of the interaction 
strength. The FLEX approximation displays the correct qualitative behavior of developing a peak, 
but is off by an order of magnitude in the peak position and height. The simplified FLEX* yields 
a quantitatively more accurate approximation, but is poorer in the limit of weak coupling because 
it does not include all of the third-order contributions properly. 
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